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UNSYMOTRICAL  BUCKLING  OF  THIN  SHALLOW  SPHERiaL  SHELLS* 

by 

Nal'Chlan  Huang 


INTRODUCTION 

This  papar  la  concarnad  with  tha  thaoratical  study  of  buckling  of  claapad 

shallow  spharlcal  ahalls  undar  unlfom  axtarnal  prasaura  (Figure  1).  For 

sufficiently  large  daflactlon,  daforaatlona  of  such  shells  are  not  proportional 

to  the  applied  pressure.  The  shell  deforms  axlsyonatrlcally  under  sufficiently 

low  pressure  and  If  wo  asstsae  that  the  deflection  naaalns  axlsynssetrlcal,  the 

pressure-deflection  relation  nay  be  represented  by  a  curve  such  as  OAB  In 

Figure  2.  At  the  suxlnus  pressure  ,  the  shell  tends  to  Jump  from  A  to 

B  .  So  q  Is  the  buckling  pressure  for  axlsyneetrlcal  snapping.  Tha  problem 
cr 

of  axlsymmetrlcal  snapping  has  been  solved  by  different  ntsierlcal  methods  In 
References  (1-4)  and  tha  results  agree  with  each  other  and  are  represented  by 
the  curve  In  Figure  3.  As  shown  In  Figure  3,  the  buckling  pressures  obtained 
In  such  a  manner  are  too  high  as  compared  with  experimental  results  obtained  in 
References  (5-6).  Initial  Imperfections  of  the  shell  and  unsynmetrlcal  buckling 
are  presumed  to  be  the  sources  of  this  discrepancy  between  axlsynmtetrlcal  buckling 
theory  and  experiment.  UnsyMetrlcal  deformation  could  conceivably  start  to 
develop  at  point  C  In  Figure  2.  After  unsynetricel  deformation  has  been 
superimposed,  the  pressure-deflection  curve  might  be  represented  by  a  branch 
CD  or  CE  at  C  as  shown.  Tha  pressure  at  C  Is  defined  as  the  critical 
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prtfture  for  uaaTaMtrical  buckling.  Grlgolyuk^^^  applied  'the  Galerkln  method 
to  aolve  the  problem  of  uaa3nBmetrlcel  buckling,  but  no  niaierlcal  raeulte  were 
given  In  hie  work.  GJelavlk  end  Bodner'^  end  Farmerter  and  Fung'  '  worked  on 
thla  problem  by  the  energy  method  end  numerical  calculations  were  carried  out 
only  for  a  particular  unsymmetrlcal  buckling  mode.  Welnltschke^^^^  used  a  power 
series  method  to  solve  this  problem  and  obtained  extensive  numerical  results, 
but,  as  will  be  shown,  the  buckling  pressures  obtained  In  his  work  are  In  serious 
disagreement  with  the  results  of  this  paper. 

BASIC  EQUATIONS  OF  SHALLOW  SPHERICAL  SHELLS 

A  shell  Is  called  "thin"  If  the  ratio  of  Its  thickness  to  the  radius  of 
curvature  of  Its  middle  surface  Is  much  leas  than  unity;  and  a  spherical  shell 
is  called  "shallow"  if  the  ratio  of  its  rise  at  the  center  to  the  base  diameter 
Is  less  than,  say,  1/8  .  The  middle  surTace  of  a  shallow  spherical  shall  can 
be  represented  by  the  paraboloid 

s  -  hCi-(;)^]  (1) 

where  H  is  the  rise  of  the  middle  surface  at  the  center  and  a  is  the  base 
radius,  as  shown  In  Figure  1.  The  radius  of  curvature  of  the  shell  Is 

*  "  2H 

The  stress  resultants  and  bending  moments  per  unit  length  of  the  shell  are 
shown  In  Figure  4,  where  N^  ,  N^  and  N^g  are  membrane  forces  respectively; 


Ntaibcrs  in  the  raised  parenthesis  refer  to  the  references  at  the  end  of 
this  report. 
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Q  and  Q.  ara  tranaverae  aheara  and  M  ^  M  and  M  are  Mrldlonal, 
r  0  XT  0  jr0 

clrcimferantlal  and  twiatlng  momenta  raapectlvely.  Let  U  ,  V  and  W  be  the 
horizontal  radial^  horizontal  tangential  and  vertical  componenta  of  dlaplace- 
ment  reapectively  (Fig.  4)  and  let  q  be  the  external  preaaure  which  can  be 
conaldered  to  be  vertical. 

In  the  following  equatlona,  we  uae  the  notatlona 

()•■!;()  and  (*)  ■  ) 


Bqulllbrltmi  of  momenta  requlrea 

*  "re  ■  "s  ■ 

<'V  *  S  *  "r8  ■  ">9  ■  “ 

Equilibrium  of  atreaa  reaultanta  requlrea 

(r»,)'  -  0  (5) 

Wt,)'  -0  «) 


+  C7  +  0,]  +  ti|  -  0 

Note  that  nonlinearity  has  been  Introduced  In  Eq.  (7)  by  considering  the 
Influence  of  W  . 

The  stress  strain  relations  are 


^re 


r?(N.-vN  ) 
tE'  S  r' 


2(l+v)  K 
tE  re 


(7) 


(8) 

(9) 

(10) 


M  -  D(h  +v  K.) 
r  r  0 


(11) 
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Mj  -  D(Kg+v  K^) 


“re  “  ^e 


(12) 

(13) 


iA«r«  I  and  v  ara  tha  alaaclc  modulua  and  Poisaon'a  ratio  raspactlvaly;  t 

3 


■Canda  for  tha  thlcknaaa  of  tha  aha 11  and  D  > 


Kt' 


12(l-v'') 

Tha  attain  diaplacamant  ralationa  ara  alao  nonllnaar 

-  O'  -  a'W  +  j(W')^ 

e.  -  -  U  +  -  V  -  -ir  iW  +  i(-  W)^ 

e  r  r  ^2  2'r  ' 

Y  .  1  6 . 1  V  +  V  -  -  iW  -  -  a'M  +  -  iw' 

're  r  r  r  r  r 

K  -  -W" 
r 

Equatlona  (5)  and  (6)  can  ba  aatiaflad  by  aattlng 

r 


“e  - 


N  .  -  -(-  F)’ 
r9  r 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 
(22) 


whara  F  la  a  atraaa  function. 

Ellalnatlng  tha  tranavaraa  ahaara  Q  and  Q.  froa  Bqa.  (3),  (4)  and  (7), 

r  0 

we  have 


[(rM^)'  +  -  Mj]'  +  j  C(rM^g)’  +  rN^(W-*)"  +  Ng[^(W-a)  +  (W-a)'] 


+  2Nj.g[(W-a)'  -  ^(W-z)]  +  q^.  •  0 


(23) 
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Thtn  the  aquation  of  oqulllbrlua  obtalnad  by  direct  aubetitutlon  is 
DV^V^W  -  (-  F'  +K  P)(W-s)''  +  ^  r'i^z)  +  -  F"(M-r)' 


+  2(^  r  -  ^  F)(j  i'  -  ju'  i  +  ^V)  +  q  (24) 


where  )  w  (  )"  +i(  )•  +  ^('')  . 

r 


(25) 


(26) 


Using  Eqs.  (8) -(10)  and  (20) -(22),  we  can  show  that 

■  r  ■  7Wm>'  *  ■  E 

Substituting  Iqs.  (14)-(16)  Into  the  left  hand  side,  we  get  another 
fundaieental  equation  Involving  F  and  cu 

A  -  Kt  |-2(^  ^  *)<r  "'*  ^  +  *’’^r  +  "^  *)  +  W"(^  z’  +  ^  t) 

+  (^  W*  -  W)^  •  M"(^  W  +  -^  if)} 

which  Is  the  coeipatlblllty  equation  of  shallow  spherical  shells. 

Substituting  Eq.  (1)  Into  Eqs.  (24)  and  (26),  we  get  Marguere's  nonlinear 
differential  aquatlons^^^^  In  polar  coordinates  for  shallow  spherical  shells 
under  unlfone  pressure, 

DvS  -  ^  V^F  +  (^  F'  +  ^  'r)\r'  +  (^  W  +  “  W)F"  -  2(^  i'  -  F)(^  W'  -  W)  +  q 

r  r  r  r 

(27) 

^F  -  Etc-  ^  V^W  +  (J  W*  -  fr)^  -  (-i  W  +  H)M"] 


(28) 


r  r 

Let  us  Introduce  the  following  nondlwenslonal  quantities: 

X  -  2C3(l-v^)]^^^(H/t)^^^ 

x-^r 

a 


6- 


w  —  W 
*  2H 


T  -  -=-5-  F 
4EH  t 


wh«r«  q 


32IH^t 


if  th«  clattical  buckling  prnasurn  of  ■  coaplote  apborlcal 


aholl  of  tha  aaaa  radlua  of  curvatura  and  thicknaaa.  Than  tha  nondlBMnalonal 
forma  of  Eqa.  (27)  and  (28)  ara 

-  V^T  +  (^  F'  +  ^  »)w"  +  (^  w'  +  \  w)'r"  -  2(^  r  -  ^  i')(^  *•  -  ~  w)  +  4p 


X 

2  2 

whara  ()'■!;:()  ‘“d  "(^  (  )  *  (^  ^  ^  ' 

bx^  *  b6^ 

The  adga  of  cha  aha II  la  completely  clamped;  therefore^  on  tha  boundary 


r  <■  a  ,  wa  have 


W  -  0 


W'  -  0 


U  -  0 


V  -  0 


From  tha  laat  two  boundary  condltlona,  wo  can  ahow  by  tha  atraln  dlaplaca- 


mant  ralatlona  that 


«e  -  0 


on  tha  boundary  r  -  a 


7' 


Expr««aad  In  tana  of  tho  atraaa  function,  thaaa  oquatlona  bocoM 

F"  -  -  F'  -  F  -  0  (37) 

a 

a(F”  -  F'  -  F)'  -  -  F'  -  ^  F  +  vF”  +  2(l+v)(^  F) '  -  0  (38) 

*  a 

on  r  ■  a  . 

Tha  nondlaanatonal  foma  of  Eqa.  (31),  (32),  (33)  and  (34)  ara 

w  -  0  (39) 

w*  -  0  (40) 

f"  -  f  f  V  -  0  (41) 

X(r  -  ^ f  -  if)'  -  ^  Y  +  vr  +  2(l+v)(^  T)'  -  0  (42) 

X  X 

on  X  ■  X  . 

In  addition,  at  tha  canter  of  tha  ahall,  w  and  T  Buat  fulfil  tha  raquira- 
ment  of  finitt;  raaaea. 


GOVERNIHG  EQUATIONS  FOR  BUCKLING  OF  CLAMPED  SHALLOW  SPHERICAL  SHELLS 

Aa  mentioned  before,  the  ahall  neceaaarily  deforma  axlaymmatrlcally  only 
under  aufficiently  low  preaaura.  The  preaaure  might  reach  a  critical  limit 
(point  C  in  Fig.  2)  auch  that  tho  ahall  bifurcataa  from  an  axiaymmatrlcal 
deformation  path  to  one  of  unaynmiatrical  deformation. 

Juat  after  bifurcation,  the  functiona  w  end  T  can  be  written  aa 


a3*(x)  +  a)(x,8) 

(43) 

♦*(x)  +  i|r(x,e) 

(44) 

where  (u*(x)  and  <|r*(x)  ara  tho  nondimanaiondl  vertical  deflection  and  the 


8- 


nondlmenslonal  atreit  function  Juat  before  buckling  and  hence  are  axlayametrical; 
<n(x;0)  and  ^(xj9)  are  due  to  unaymetrlcal  buckling  and  theae  are  considered 
to  be  inflnltealmally  amall  at  the  beginning  of  unaynBetrlcal  buckling.  From 
Bqa.  (29)  and  (30),  we  have 

^0^  -  ^  ^*"05*'  +  4p  (45) 

.  i  (46) 

Put  e*  ■  -oJ*'  and  f*  ■  ♦*'  ,  Eqa.  (45)  and  (46)  can  be  written,  after 
one  integration,  aa 

(x6*')  -  +  xf*  -  -2px^  +  0***  (47) 

(x**')'  -  ^  -  x0*  -  -  I  0*^  (48) 

Alao,  the  boundary  eondltlona  can  be  reduced  to 

0*(X)  -  0  (49) 

Xf*’ (X)  -  v?*(X)  -  0  (50) 

Theae  are  the  governing  equations  of  the  syanetrlcal  problea  which  has  been 
solved  numerically  In  References  1-4. 

Substituting  Eqa.  (43)  and  (44)  into  Eqa.  (29)  and  (30)  and  using  Eqa.  (45) 
and  (46) ,  we  obtain 

-  (-  t'  +  V)9*'  +  -  -  -  y"0*'  +  (-  0)'  +  -^  c6)0*'  (51) 

X  X 

A  ■  ^x 

X  * 

whore  only  the  linear  terms  In  cc  and  y  have  been  retained.  These  are  two 
linear  differential  equations  for  the  functions  co  and  y  .  Hence  our  problem 
has  been  reduced  to  an  eigenvalue  problem  where  the  eigenvalue  p  is  implicit 
and  involved  in  the  solution  0*  and  of  the  nonlinear  problem  as  defined 
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by  Eqs.  (47) -(50). 
Let 


m(x 


>«>  -  z  % 


(x)coa  n9 


n“0 


(53) 


♦  (x^e)  -  ^  \|r^(x)coe  ne  (54) 

n-0 

2  2 

n  .  2  X  ux  2  n  n  n 

dx  X 

Substituting  Bqs.  (53)  and  (54)  Into  (51)  and  (52)^  ws  have,  for  each  n  , 

2  2 

L^OJ  -  L  *  -  (^  ♦ '  -  %  ♦_)e*'  +  ^  tu"**  -  ^  t''9*  +  ‘“n  '  0)  )♦*'  (55) 

n  n  n  n  x  n  2  n  3C  n  x  n  x  n  2  n 

X  X 


L^i  -  -L  CO  +  (-  (O'  -  CO  )0*'  +  -  co"e* 
n  n  n  n  x  n  2  n  x  n 

X 


(56) 


The  boundary  condition*  can  be  derived  from  Eq*.  (39)-(42).  The**  ere 


co^(X)  -  0 


co'(X)  -  0 
n 


(57) 

(58) 


2 

♦  "(X)  -  Ti’M  +  "  0 

n  ^  n  ,  ^2  n 


(59) 


Xr'(X)  -  7  [1  -  V  +  (2+v)n2>'(X)  ♦„(X)  -  0 

n  .  X  n  1  n 


(60) 


Equation*  (55) -(60)  are  th*  governing  differential  equation*  and  the 
boundary  condition*  of  the  problem  of  uneynwtrlcal  buckling  of  clamped  ahallow 
apherical  ahell*.  Exlatance  of  a  aolutlon  for  any  Integer  n  at  a  preaaure  lea* 
than  the  critical  preaaure  of  axlayaaaetrlcal  snapping  Inpllea  exlatence  of  a 
critical  preaaure  for  unaynmetrlcal  buckling. 
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SOME  PRELIMINARY  STUDIES 


Modified  Probleon 

Before  the  mein  problem  wee  treated.  It  wee  considered  advisable  to  study 
some  problems  of  unsymnetrlcal  buckling  of  shells  with  modified  boundary 
conditions  for  which  the  solutions  can  be  obtained  analytically.  The  purpose  of 
this  study  Is  to  find  whether  unsymnetrlcal  buckling  rather  than  axlsymmetrical 
snapping  would  control  In  these  modified  problems.  The  answers  for  the  critical 
pressure  for  unsymmetrlcal  buckling  obtained  here  can  be  also  used  to  check  the 
accuracy  of  the  numerical  procedures  for  the  main  problem. 

In  the  first  modified  problem  the  edge  of  the  shell  Is  supported  by  rollers 
which  can  slide  along  a  conical  wall  without  changing  the  slope  along  the  edge 
of  the  shell  as  shown  In  Figure  5.  In  the  second  modified  problem  the  edge 
condition  Is  the  same  as  that  In  the  first  modified  problem  except  that  the  edge 
Is  Imagined  to  be  suddenly  clamped  Just  before  buckling  occurs.  In  both  modified 
problems  the  shell  contracts  uniformly  before  buckling  and  therefore,  0*  ~  0 
and  i*  ■  -2px  . 

The  critical  pressures  for  these  modified  problems  are  calculated  for 
n  ■  0,  1,  2  and  3  In  Appendix  A  and  are  plotted  against  \  In  Figures  5  and  6, 
where  n  la  the  number  of  waves  along  the  circumferential  direction  In  the 
buckling  mode.  In  both  modified  problems  It  is  found  that  unsymmetrlcal  buckling 
(n^O)  dominates  for  some  values  of  X  .  These  calculated  critical  pressures 
are  higher  than  the  classical  buckling  pressure  of  a  complete  spherical  shell. 

In  the  first  modified  problem  the  calculated  critical  pressures  meet  the 
classical  buckling  pressure  periodically  at  certain  values  of  X  and  In  the 
second  modified  problem  the  calculated  critical  pressure  approaches  the  classical 


-11- 


buckllng  pressure  asympCoticelly  when  X  approach  Infinity.  In  both  problems 
the  results  are  not  realistic  because  of  initial  imperfections. 

Approximation  by  Variational  Principle 

A  variational  principle  is  developed  in  Appendix  B  for  a  preliminary  study 
of  the  critical  pressure  of  unsyssaetrlcal  buckling  of  clamped  shallow  spherical 
shells.  The  boundary  conditions  (57)  and  (58)  of  the  main  problem  were  specified 
in  the  application  of  this  principle.  The  differential  equations  (55)  and  (56) 
are  the  Euler 'a  equations  and  the  boundary  conditions  (59)  and  (60)  are  the 
natural  boundary  conditions  of  this  variational  principle.  A  Rayleigh-Ritz 
method  is  used  in  conjunction  with  the  variational  principle  in  Appendix  B  and 
numerical  results  have  been  obtained  for  X  ■■  6  as  shown  in  Table  1.  These 
critical  pressures  for  unsymmetrical  buckling  are  lower  than  the  critical  pressure 
for  axisymmetrical  snapping^  and  this  certainly  tends  to  confirm  the  suspicion 
that  unsymmetrical  buckling  does  occur  for  clamped  shallow  spherical  shells. 
However,  the  numerical  results  obtained  by  this  variational  principle  are  not 
necessarily  either  upper  bounds  or.  lower  bounds  to  the  exact  critical  pressure 
for  unsymmetrical  buckling.  . 

NUMERICAL  SOLUTION 

Difference  Equations 

The  most  popular  practical  method  to  solve  differential  equations  is  to 
reduce  them  to  difference  equations  by  confining  the  range  of  the  independent 
variables  to  a  network  of  mesh  points.  In  this  problem  we  define 
x"jh  J«0,1,...,N»  XAi 

where  h  is  the  mesh  size. 


wrlctsn  88 
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Vj-1  "  °  J  -  1,  2,  ...  N  (71) 

where  one  fictitious  station  J  -  IH-1  has  been  added  off  the  edge  of  the  shell 
for  convenience.  Equations  (65)- (68)  can  be  written  as 

■  Gy„.i  ■  0  (72) 

where 


-14 


(82) 


1 1 
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From  Sq.  (79) 

Ql  -  (83) 

which  also  can  ba  obtalnad  from  Eq.  (82)  by  setting 

Qq  -  0  (84) 

Equation  (82)  la  a  racurranca  formula  by  which  all  Q's  can  ba  avaluatad. 

Substituting  Eq.  (80)  Into  Eq.  (72),  wa  gat 

Sy^l-0  (85) 

whara 

S  -  G(I-Qj,.^Qj,)  +  KQ^  (86) 

and  I  Is  a  4x4  unit  siatrix.  For  a  nontrivial  solution  G  , 

|S|  -  0  (87) 

which  Is  the  characteristic  aquation  for  determination  of  critical  pressures  of 
unsymmetrlcal  buckling. 

The  case  n  ■  1  Is  a  special  case  where  Eqs.  (55)-(60)  can  ba  satisfied  for 
all  p  by  an  exact  solution  -  0  and  *  x  .  Since  the  difference 

equations  are  exact  fot  linear  functions  Eqs.  (71)  and  (72)  can  bo  satisfied  by 

0 

J 
0 
0 

Equation  (85)  can  be  written  as 

o' 

[s]  -  0 

0 


0 
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for  n  ■■  1  and  to  all  elements  in  the  second  column  of  matrix  S  vanish.  For 
n  -  1  ,  Eq.  (56)  Is  reduced  to 


4-  [x\"'  -  3\(r'  +  —  i|f  +  x^co'  -  xo)  -  (xa)'-a))fl]  ■  0 

dx  X 

Integrating  from  x  0  to  x  *  X  and  using  Eqs.  (57)  and  (58)^  we  can  obtain 
Eq.  (60).  ThuSj  for  n  >  1  ^  Eq.  (60)  Is  dependent  and  can  be  omitted.  We  have 


[S] 


■“Wf 

“wfl 


0 


iriiere  S  la  the  resulting  matrix  obtained  by  striking  out  the  second  column  and 
the  fourth  row  of  the  matrix  S  .  For  a  nontrivial  solution, 


|S|  -  0 


(88) 


which  Is  the  characteristic  equation  for  determination  of  critical  pressures  for 
unsymmetrlcal  buckling  for  n  •  1  . 

Numerical  Procedure  and  Results 

The  numerical  procedure  of  this  problem  consists  of  two  parts.  The  first 
part  Is  to  calculate  6*  ,  and  their  derivatives  for  the  axlsymmetrlcal 
deformation  before  unsymmetrlcal  buckling.  We  used  the  Iterative  method  given 
In  Reference  (1)  .  The  Poisson's  ratio  v  was  taken  equal  to  1/3  In  all 
calculations  and  the  Interval  tOi  used  was  0.25.  The  first  derivatives  of  6* 
and  f*  were  calculated  by  Eq.  (69).  The  values  of  Kelvin  functions  and  their 
first  derivatives  used  were  found  In  Reference  (12). 

In  the  second  part  of  the  calculation,  the  elements  of  matrices  ,  Bj  , 


The  criterion  of  convergence  of  the  Iterative  procedure  used  in  this 
computation  Is  that  the  final  answer  for  average  deflection  agrees  to  four 
significant  figures  with  the  mean  of  the  results  of  the  five  previous 
Iterations , 
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Cj  ,  G  and  K  wera  calculated  first  for  different  values  of  n  and  ,  S 
and  |s|  could  be  obtained.  In  general^  |s|  does  not  vanish  If  the  trial  value 
of  p  Is  not  the  buckling  pressure  values  of  |s|  were  plotted 

against  p  .  The  critical  pressure  Is  the  lowest  root  of  |s|  ■■  0  .  When  n  ■  1 
S  must  be  calculated  Instead  of  S  .  This  procedure  continued  for  different 
values  of  n  until  axlsymetrlcal  snapping  occurred  before  unsymmetrlcal  buckling 
appeared . 

Two  nesh  sizes ^  h  ■  0.25  and  h  0.125  ,  were  used  In  the  calculation  for 
X  ■  6  .  Since  the  nesh  size  0.25  was  used  In  the  first  part  of  the  calculation. 

In  the  computation  for  h  >  0.125,  the  values  of  6*  ,  i*  and  their  first 
derivatives  at  the  Intermediate  stations  were  obtained  by  linear  Interpolation. 

The  difference  of  buckling  pressures  calculated  by  these  two  mesh  sizes  was  about 
0.1  %  .  Hence  we  use  h  «  0.25  In  our  calculation  . 

The  validity  of  this  numerical  procedure  was  checked  by  comparing  solutions 

of  the  second  modified  problem  In  Appendix  A  (Fig.  6)  with  solutions  obtained  by 

this  numerical  method .  Computations  were  carried  out  for  the  cases  X  ■■  6  and 

X  ■  10  and  the  maximum  error  In  p  was  0.4  %  . 

cr 

All  numerical  calculations  were  made  on  the  IBM  7090  digital  computer  and 
the  prograimnlng  was  written  In  the  FORTKAH  language.  The  critical  pressures  for 
different  values  of  X  and  n  are  given  In  Table  2  and  also  plotted  In  Figure  8. 
The  lowest  value  of  these  critical  pressures  at  a  given  X  Is  the  governing 
pressure  for  unsymetrlcal  buckling.  These  governing  pressures  are  shown  by  a 
heavy  line  In  Figure  7..  Buckling  modes  were  evaluated  for  these  governing 


This  choice  of  mesh  size  provides  a  reasonable  number  of  stations  within 
the  boundary  layer  of  the  buckling  mode. 
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pressures  and  are  ahown  In  Figure  8. 

According  Co  the  reaulta  of  Che  calculaClon^unsynoaeCrlcal  buckling  (n  ^  0) 
does  noc  occur  for  X  <  5.5  .  UnaymneCrlcal  buckling  aCarCa  Co  appear  when  \ 
la  allghcly  greaCer  chan  5.5.  As  X  keepa  Increaalng  Che  buckling  mode  ahowa 
more  and  more  wavea  along  Che  clrcionferenClal  dlrecClon  and  alao  ahowa  a  dlaClncC 
boundary  layer  near  Che  edge  of  che  ahell  along  Che  radial  dlrecClon  when  X  la 
high.  An  aaympcoClc  analyala  haa  been  done  In  Appendix  C  where  an  aaympCoClc 
value  of  Che  crlClcal  preaaure  for  unaymmeCrlcal  buckling  la  found  Co  be  0.864 
when  X  approachea  InflnlCy,  and  Che  raClo  n/X  la  found  Co  approach  0.817. 


DISCUSSION 

The  ahell  may  anap-Chrough  under  Che  crlClcal  preaaure  for  unaymmeCrlcal 
buckling  If  Che  cangenC  of  Che  branch  on  Che  preaaure  deflecClon  curve  haa  a 
negaClve  alope  aC  chaC  crlClcal  preaaure  (aa  ahown  by  Che  branch  CD  In  Fig.  2). 
On  Che  ocher  handj  if  ChaC  branch  haa  a  pbalclve  alope  (aa  CE),  deformaClon  of 
Che  shell  changes  from  che  axlsymmeCrlcal  Cype  Co  Che  unsymmeCrlcal  Cype  suddenly 
under  Che  crlClcal  pressure  buC  no  snapping  appears.  The  deCermlnaClon  of  Che 
branch  of  pressure-deflecClon  curve  Involves  Che  analysis  of  che  poac-buckllng 
behavior  of  che  shell. 

In  Figure  9  Che  resulCs  of  previous  aCCempCs  Co  calculaCe  crlClcal  pressures 

for  buckling  of  shells  are  ploCCed  for  Che  purpose  of  comparison  wlch  Che  presenC 

(8) 

resulCs.  Gjelsvlk  and  Bodner'  '  and  ParmerCer  and  Fung'  '  calculaCed  Che 
crlClcal  pressures  based  on  an  approxlmaCe  soluclon  for  Che  cases  n  ■  0  and 
n  ■  1  .  The  Juncclons  of  curves  for  n  0  and  curves  for  n  ~  1  In  Figure  9 
are  represenCed  by  clck  marks.  Helnlcschke^^^^  obCalned  crlClcal  pressures  for 
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■n  axtenslvc  rang*  of  X  and  n  by  a  method  which  waa  claimed  to  be  accurate, 

but  the  buckling  presaurea  obtained  In  hla  work  are  In  aerloua  dlaagreement  with 

the  preaent  reaulta.  It  la  noted  that,  for  n  >  1  the  preaent  reault  la  much 

* 

cloaer  to  the  reaulta  of  Referencea  (8)  and  (9)  than  Welnltachke  a  reault  .  In 
Referencea  (7),  (8)  and  (9)  the  curve*  repreaentlng  critical  preaaurea  for  n  •  0 
are  ahown  to  be  tangent  to  thoae  for  n  -  1  at  their  Junction*.  However,  thla 
tangency  doea  not  appear  In  the  curve  of  preaent  work.  The  preaent  theoretical 
buckling  preaaurea  are  atlll  higher  than  the  experimental  reaulta.  The  effect* 
of  Initial  unaymmetrlcal  geometrical  Imperfection*  and  variation  of  ahell  thlckneaa 
are  preaumed  to  be  the  aource  of  thla  dlacrapancy,  but  the  analyala  of  aucb 
problem  la  vary  comi^llcated. 


More  recently.  In  a  private  communication,  G.  A.  Thuraton  atated  that  he 
found  a  lower  bound  to  the  critical  preaaure  p^^  ■  Q.753  for  n  ■  4 

and  X  ■  8  ,  which  la  cloae  to  the  preaent  reault  p^^  ■  0.766  (Table  2). 
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APPENDIX  A 

Analysis  of  Shells  with  Modified  Boundary  Conditions 
As  mentioned  before^  In  both  modified  problems 

e*  -  0 

i*  ■  -2px 

Substituting  them  Into  Eqs.  (55)  and  (56)  one  obtains 

2 

L  (O  •  L  *  -  2pL  a> 

n  n  n  n  n  n 


L>  -  -L  uj 
n  n  n  n 

The  general  solutions  are  expressed  In  terms  of  Bessel's  functions 


(Al) 

(A2) 

(A3) 

(A4) 


05  -  A  J  (^ix)  +  B  J  (-)  +  C  x"  +  E  Y  (^x)  +  F  Y  (-)  +  G  x'" 
n  nn  nn^t  n  nn  nnn  n 

A  ,  E  5 

n^n  nnu  n  ^n  nnu  n 


where 


p  -  I  (n^+  z^) 


(A5) 

(A6) 

(A7) 


All  components  of  stress  resultant  and  moment  are  bounded  at  the  center  of  the 
shelly  therefore, 

05  -  A  J  (px)  +  B  J  (-)  +  C  x" 

XI  nn  nn^i  n 

♦  -  J  (nx)  +  B  (f;)  +  D  x" 

n  n  n  n  4  n 

The  boundary  conditions  of  the  first  modified  problem  are  conveniently 

obtained  from  the  principle  of  virtual  work 
2x  a 

drdS  -  0 

0  0 


(A8) 

(A9) 


(AlO) 


-21- 


Uslng  the  atreln  diepleceMent  relatione  we  get 
2*  *  ^ 

0  0 

+  C(rN^)  '+N^g-N0]6U  +  [(rN^g) 

2k 

+  /  /C(rM^)'+2M^g-Mg]6C  +  ((rM^) '+2M^g-Mg+2HN^)+aN^]6?  -  aM^6w' 

0  ^ 

+  «N  *6vl  de  -  0  (All) 

®  J  r-a 

where  C  en<l  i  ere  the  dlaplacaeianta  along  the  wall  and  perpendicular  to  the 
wall  reapectively,  therefore, 

^  M  u  .  w  (A12) 

(A13) 

(A14) 

(A15) 
(A16) 
(A17) 

(A18) 
(A19) 
(A20) 


+  (l+v)(n^-l)*  +0)^-0 

n  n  n  n  n  n 


?  -  U  +  —  W 

’  a 

By  aubatltutlon,  we  get  the  boundary  condltiona 

(rM  )'  +  2M  -  -  M-  -  0 
r'  r0  6 

W  -  0 

?  -  0 

N  .  -  0 
r0 

In  tenaa  of  the  nondlaienalonal  quantltlea,  theae  are 

co'  -  0 
n 

2 

n  n  ,  4  n 


(A21) 
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Substituting  Eqa.  (A8)  and  (A9)  into  Eqs.  (A18)-(A21)  and  allainatlng  tha 
coefficient  of  ,  we  have  the  following  characteristic  equation  froa  the 
resulting  equations 


^  J'(^) 
H  nV 


^l\V(^iX)  +  (l-v)nV(MX)  -  ^  j;(J)  +  (l-v)n^J„(^)  (l-v)n^  |-  0  (A22) 


ix  -  ±2 

X 


^  -  n  \  Jjb  ^ 


X  n  V 


2  “'aV' 


(l+v)X" 


The  case  n  ■  1  is  also  a  special  case  in  the  modified  problems.  However, 

Eq.  (A22)  still  holds  for  n  -  1  . 

The  buckling  pressures  are  plotted  against  \  for  different  values  of  n 
in  Figure  5.  The  case  n  >  0  represents  an  axlsynmetrical  buckling. 

The  boundary  conditions  of  the  second  modified  problem  are  the  same  as  those 
of  the  main  problem.  Substituting  Eqs.  (A8)  and  (A9)  into  Eqs.  (57)-(60)  and 


eliminating  the  coefficient  of  ,  we  get  the  following  characteristic  equation 


J  (-) 

n'ji 


HXJ^(HX) 


^  J'(-) 

H  n  n 


S:  j  (lix)  f  j'(i)  -  j  A  -  . 

HA  n  2n^X^  X  n  )i  n  |i  1+v 

This  equation  holds  for  the  special  case  n  -  1  but  is  not  applicable  for  n  ■  0 
idiere  we  can  show  by  using  Eqs.  (55)  and  (56)  that  Eqs.  (57)-(60)  are  linearly 
dependent.  Omitting  the  last  boundary  condition  equation  (60),  we  can  get  the 
characteristic  equation  for  n  >  0  ,  which  is 

t  +  (l+v)(n2.  -  0  (A24) 


The  p  -  X  curves  for  the  second  modified  problem  are  shown  in  Figure  6. 
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AFPBNDIX  B 

Variational  Principle  for  Buckling  of  Clamped  Shallow  Spherical  Shells 

A  variational  principle  for  the  problem  of  clamped  shallow  spherical  shells 

(13) 

was  first  considered  by  Weinitschke  for  the  problem  of  axisymmetrical 

(14) 

deformation  of  shells.  From  a  modified  Reissner  variational  principle  ,  we 
can  show  that  Eqs.  (45),  (46),  (49)  and  (50)  are  equivalent  to  the  following 
variational  principle: 

6  / {x  j(^V*)^-4pa^}  xdxf  j(l+v)6Ct*'.(X)]^  -  0 


where  a:^(X)  ■  a3*'(X)  "0  ia  specified.  Furthermore,  Bqs.  (55) -(60)  era 
equivalent  to  the  following  variational  principle: 

6  t*'(<u')^+ 7  t*"  )(xa>’-kd*'a)')-n^(- ♦  )  ' (o)  +  ' 

J  1 2x  ^  '  n*'  2  ^  2'  n  2  n  n  ^x  n  x  ' 

^  X  X 


where  w  (X)  ~  oo'  (X)  *0  is  specified  and  variations  are  taken  with  respect  to 
n  n 

cu  and  i|r 

In  order  to  apply  these  variational  principles  we  try 

CJ*  .  (x^-X^)(Ax^+BX^)/x'^  (B3) 

♦*  -  Cx^  +  DX^x^  (B4) 

0)  -  Ex"(x^-X^)^  (B5) 

n 

n  4  2  2  4 

-  X  (Fx  -HJX  X  +HX  )  (B6) 


Substituting  these  expressions  into  Eq.  (Bl)  and  taking  variations  with 
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reapect  to  A  ^  B  ,  C  and  D  ,  wa  obtain  four  nonlinear  equatlona  for  the 
unknowns  A  ,  B  ^  C  and  D  ;  their  values  can  be  calculated  numerically  for 
any  assigned  p  .  Substituting  Eqs.  (B3)-(B6)  into  Eq.  (B2)  and  taking 
variations  with  respect  to  E  ,  F  ,  6  and  H  ,  we  obtain  four  homogeneous 
linear  equations  and  hence  the  characteristic  equation  from  which  the  buckling 
pressure  can  be  determined  for  any  n  .  The  numerical  calculation  has  bean  done 
for  X  "  6  and  the  buckling  pressures  are  shown  in  Table  1. 
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APPENDIX  C 


Aayaptotlc  Solution  for  Buckling  Pressure  for 
Shallow  Spherical  Shells  with  Large  X 
It  Is  Interesting  to  consider  the  asymptotic  behavior  of  clamped  shallow 
spherical  shells  when  X  approaches  Infinity.  A  boundary  layer  Is  found  near 
the  edge  of  the  shell  i^en  the  shell  deforms  axlsysmietrlcally  before  buckling. 

In  the  region  outside  the  boundary  layer  the  shell  deforms ^essentially  by  a  rigid 
body  downward  displacement^  hence  all  components  of  stress  resultant  and  moment, 
except  N^  ,  vanish. 

Ut 

i  -  {*  +  2px  (Cl) 

Equations  (47) -(50)  csn  be  written  as 

I*"  ^  (“) 


(C3, 

X 

§*(X)  -  0  (C4) 

f'(X)  -  ^  l(X)  -  2(l-v)p  (C5) 

In  the  region  of  boundary  layer  where  x  approaches  Infinity,  the  following 
equations  can  be  used: 


9*"  +  i  -  -2p9*  (C6) 

i"  -  0*  -  0  (C7) 

e*(X)  -  0  (C8) 

i'(X)  -  2(l-v)p  (C9) 
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and 


Alto 

03(0)  ■  a3''(0)  ■  \|f(0)  ■  t"(0)  ■  0  (C21) 

By  the  same  method  as  given  in  the  main  problem  the  asymptotic  value  of 
the  critical  pressures  can  be  evaluated  numerically.  In  the' calculation^  6^ 
was  chosen  equal  to  40.  The  critical  pressures  p^^  are  plotted  against  o 
in  Figure  10.  The  minimum  of  this  curve  determines  the  required  asymptotic  value 
of  the  critical  pressure  for  X  equal  to  infinity  which  is  found  to  be 
P^^  ~  0.864  when  the  ratio  n/X  approaches  0.817.  The  buckling  mode  is  also 
calculated  and  is  shown  in  Figure  11  from  which  the  actual  boundary  layer 
thickness  is  found  to  be  12  approximately. 
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FIG.  2  PRESSURE -DEFLECTION  CURVE  OF  CLAMPED 
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AND  EXPERIMENTAL  RESULTS 


FIG.  4  STRESS  RESULTANTS,  MOMENTS  AND 
DISPLACEMENTS  IN  SHELLS 


SECOND  MODIFIED  PROBLEM 


FIG.  7  CALCULATED  BUCKLING  PRESSURES  OF  CLAMPED  SHALLOW 
SPHERICAL  SHELLS  AND  EXPERIMENTAL  RESULTS 


FIG.  8  BUCKLING  MODES  OF  CLAMPED  SHALLOW 
SPHERICAL  SHELLS 


FIG.  9  COMPARISON  OF  CALCULATED  PRESSURES  OF  UNSYMMETRICAL 
BUCKLING  WITH  PREVIOUS  RESULTS 
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Olrooto,' 

Movol  Roooorch  Lohoratory 
Woohincton  25,  D4C, 


ConMnllhf!  Offlcor 

CIIR  Broinoh  Offleo 

John  Ororor  Ubrory  Bulldlnr 

66  toot  Randolph  StrMt 

Chloo«o  U,  nilflola  (1) 

CooMivllni  Offloor 

Offloo  of  Royal  Roooanh 

lu-aneh  Offloo 

1000  Ooary  Sbraot 

9an  Pranelaee,  California  (1) 

Dlrootor 

Haval  Roaaaroh  Laboratory 
Vaihinitan  25,  D.C* 

Attni  Coda  6200  (1) 

Dlrootor 

Haval  Roaaaroh  Laboratory 
Vaahlngtcn  25,  DtOt 


Suporiatondant 

U.a.  Haval  frovlnf  Drounda  SJIm  (l) 

i>aM<ron,  Vlrflnla  (i)  Wuhlniton  25,  D.O. 


teoal  tedikanea  tait  Station  Haval  Ordnanea  loot  Statto 

S,.,.™,  o».-  Ul., 

Attni  Moehanloa  Iraneh  (1)  ^ttni  9trueturaa  Dlvlaian  (1) 
OoaaanllAa  Offloav  anl  Dlrootor  IKiporlntoadont 

EP-"  (d 

Aanapolla,  Maryland  (1) 


CoaMj^W  Oonaral 
C.9.  AlrPoroo 

Hb^l^on  25,  D«D* 

AtinilMoaroh  and  DovalcmM 


Cewndor 

Xaval  Ordoanoo  Toat  SUtion 
]^oham,  China  Late,  Callfonln 
Attni  Hiytlno  Olrlal^n  (1) 

Nr.  K4H4  XoopMn,  9aorotai7 
WoldlAC  Rooatidi  Oounall  of 
Tho  InclMorinc  foundatUn 
25  W,  J5tii  StPMt 
lov  Xcrk  16,  Ko«  Xork  (2) 

c^oHadant 
Narlna  Corpa  School* 
ftiantleo,  Vlrflnla 
Attni  Dlrootor,  Narlna  Gorpa 
Daralofnant  Oantar  (1) 

Dlrootor  of  bioUlconoo 


Offloa  of  Taehnloal  Sorvloaa 
Dapartmnt  of  CoriMroa 
Waahlntton  25,  D.O,  {3 


Offloa  of  tha  Saorotary  of  Oafamt  OoMandar 


Raaaareh  and  DavalopMant  Dlv* 
(1)  Tha  Pantafon 

Wuhlnaten  25,  D.Ct 

Attni  Taahnloal  Library  (1 


Offloa,  Sooratary  of  tha  Arvy  Chlof  of  Staff 


Tha  Pantafon 
Vaahln«ton  25,  D.O. 
Attni  inv  Library 


Dap4irtMnt  of  tha  Angr 
Vaahln|ton  25,  D.O. 

(1)  Attni  OavalopMRt  Branch 
(R  aid  D  Dlv.}  ( 


Air  Raaaaroh  aid  otvalopi.CoM 
P.O.  hm  1395 
lalilMTo  3,  Maryland 
Attni  RDWI  U) 

Dlraotor 

National  Suroau  of  atandarda 
WaahlnttOB  25,  D.O. 


U.S.  Air  Poroa  Inat.  of  Toctmol.  Kaadquartara,  0.9.  Air  Peroo 
Wrlfht-HtUraon  Air  Pcroa  W^lnrton  2^  D.O, 

Dartw.  Ohio  (i)  Brtnoh 

Attni  Chlaf,  Appl.  Naohanloi  te.  (Air  TAriota  Dlvlalon)  (1) 
U.i.  AtMlo  kMrry  Ooonlaiian  Praf.  Wiltar  I,  Danlala 
.  Waahincton  25,  0.0.  »«ho^  «£  Md 

1  . .  V  /«»  Arohttaotura 

Attni  Dlrootor  of  RaaoaMh  (2}  Unlvoralty 

Waahlneton  1,  D.O.  (1) 


Dlraotor 

Naticnal  Buroau.  of  9tandarda 
WaahWtm  >5,  D.O. 


Dlrootor 

Naticnal  Baroau  ef  atanlarda 
Waahlitcton  25,  D.O. 

Attni  Aircraft  Btruoturoa  (1) 


Chlaf  of  Staff 
Dapaitnant  of  tha  Ariy 
Waahlnfton  25,  D.C. 

Attni  Spoolal  Vaapcna  Dranoh 


Offloa,  Chlaf  of  ^Inaara 
Dapartaant  of  tha  Ingr 
Vaihln^ton  25,  D.C. 

Attni  EH0>UD  Plannlnc  Dlv., 
Civil  Worka 

Offloa,  Chlaf  of  Enflnaora 
Dapartaant  of  tha  Amy 
Waahlncton  25,  D.C. 


Br.j  tn^r  Attni  BN0*HS  Spae.  &icr. 


Offica,  Chlof  of  Ordnanea 
OopartMnt  of  tha  Arogr 
VaMilneton  25,  D.C. 

Attni  Roaaarch  and  Matarl^ 
Dranoh  (Old  R  and  D  Dlv.)  (1) 

CoMaandlng  Offle'or 
frankford  Araanal 
toidaiburc  Station 
Philadolpnla  37,  PonnaylvanU 
Attni  UWatory  Olvijion  (1) 

Dlr«tor 

Vataiwaya  BxparlMnt  Station 
P.O.  Bon  631 
Vlekaturi,  Hiialaalppl 
Attni  Library  ' 

Chlaf,  Buriau  of  Shipa 
Dapartaant  ef  tha  Navy 
VaahlACten  25,  D.C. 

Attni  Coda  371  (1) 

Chlaf,  Baroau  of  Shipa 
Dapartnant  of  tha  Navy 
Vashlneton  25,  D.C. 


Chlaf,  Buraau  of  Aaronautlca 
Dapartaant  of  tha  Navy 
Vaahintton  25,  D.C. 


Chlaf,  Buraau  of  OrdMnea 
DapartiMnt  of  tha  Navy 
Vaahlacton  25,  D.C. 


Chlaf,  Buraau  of  Ordnanea 
Dapartaant  ef  tha  Navy 
Vaahlacton  25,  D.C. 

Attni  SP 


Offlea,  Chlaf  Signal  Offlcor 
Oopartnont  ef  tha  Any 
Vadiineten  25,  D.C, 


Offloa,  Chlaf  ef  Cnclnaara 
Dapartaant  of  tha  An^ 
Vaahlncton  25,  D.C. 

Attni  CNO-NL  Library  Br., 

Ada.  Sar.  Dlv.  (1) 

Offlea,  Chlaf  of  &<Alnaara 
Dapartaant  of  tha  Amgr 
Waihlnytcn  25,  D.C. 

Attni  EN9-E8  Prot.  Conat.  Br., 
Ener.  Dlv, ,  Mil  Conat.  (1) 

CoaaiAdlnc  Offlear 

CiiCinaar  Raaaareh  Davalopaant 

Laboratory 

Port  Dalvolr,  Virginia  (1) 


Comandlng  Officer 
Jatartewn  Aroanal 
Watertown,  Kaaaachuaatta 


OtfManJaat,  U.S.  Coaat  Duard 
1300  t  Straat,  N.V. 

Waahlngtcn  25,  D.C. 

Altai  ChUf,  Taattoi  wd 
Oavaiopaani  Dlvlalon  (1} 

Director 

Langley  Aaronautloal  lab. 
Langlay  ruid,  Virginia 
Attni  Stnictnraa  Dlvlalon  (2) 


National  Solancaa  Poundatlcn 
1520  K  Straat,  k.V. 

Waahlngtcn,  D.C. 

Attni  biglnaarleg  Solaneea 

Dlvlalon  (1) 

Profaaaer  R.L*  Blepllnghoff 
Dapartaant  of  Aoronautlcal  Ebig. 
Maaaaolutaatta  IDiatltuta  ef 


U.S,  Jtorltlaa  AdalBlatratton  Ifrf'V*  ,  „  „ 

Oarnral  Adnlnletratien  Offloa  Bl«.  1512  H  Street,  K.V. 
Waahlnrtw  25,  D.C,  Vaablnfton  25,  D.C. 

AtUi  Chief.  DivlaioB  of  Prallala  Attni  teada  md  Strueluraa 
nary  Oaalgn  (D  Divlalco  (2) 


Dlraotor 

Poraat  Produeti  taboratoiy 
Madlaon,  Wlaoenaln 


Civil  Aaronautlof  Adalniatratlon 
Dapartaant  of  CoMarca 
Vaablnctm  25,  D.C. 

Attni  Chief,  Alrfram  and 


Offica  of  Ordnance  Raaaareh 

2127  IVrUa  Drive 

Duka  Station 

Ourhan,  North  Carolina 

Attni  (ilv.  of  Cng.Sclancaa  (1) 


CcnMMlng  Officer 
Squler  Signal  Laboratoiy 
Port  Konnouih,  New  Jaraay 
Attni  Coaipcnanta  and  KatarUla 


Chlaf,  Ouraau  of  Shipa 
Dapartaant  of  tha  Navy 
Vaahingten  25,  D.C, 


Chlaf,  Buraau  ef  Shipa 
Daparlnnt  of  tha  Navy 
Vaahlngton  25,  D.C. 

Attni  Coda  li23  (1) 

Chlaf,  Buraau  of  Aaronautlea 
Departaant  ef  Ua  Navy 
Uaahinglon  25,  D.C. 

Attni  AIX-2  (1) 

Chief,  Buraau  of  Aaronautlea 
Dapartaant  of  tha  Navy 
Waihlngton  25,  D.C. 

Attni  n-L2  (1) 


Chlaf,  Ouraau  of  Ordnanea 
Dapartaant  of  tha  Navy 
Vaalilnctoi  25,  D.C. 


ChUf,  Ouraau  of  Ship# 
DapartMnt  of  tha  Navy 
Waihlngton  25,  D.C. 

(1)  Attni  Coda  327 

Chlof,  Uireau  of  Shipa 
UopartMnt  of  Ua  Navy 
Wadilngton  25,  D.C. 

(1)  Attni  Code  titiO 


Profaaaer  O.f.  Carrier 
Pierce  Hall 
Harvard  Vnlvaralty 
Canbrldca  36.  KajaaehuMtta  |  I) 

Prefaater  A.C.  Crlngan 
Dapartaant  ef  Aermeutioal 
SeglAaerlng 

Putdu*  Vnlvaralty  ,  ^ 

Ufayatta,  IndUna  (1) 

Prefaaaer  L.C.  Ooudaan 
&iglnaarlng  Separiaant  Statlcn 
Vnlvaralty  ef  Rlnnaaota 
Hlnnaapolla,  Klnnaaota  (1) 


Natlmal  AcaiaiV^  of  Solenea 
21Ca  Cmtltutleat  Averai* 
Waahlngton  25,  D.C. 

Attni  taehnloal  Director,  Cai^ 
on  Shipa'  Slivotural  Daalgn  (1} 

Prefaaaer  H.H.  Dleloh 
Dapartaant  of  Civil  Zhtlnaaring 
CeluabU  Unlvoralty  ,  ^ 

New  Xerk  27,  N.T.  (1) 

rrefaaaor  Harlart  Daraaiawlea 
Northvaatam  Uolvaraliy 
Bvanaten,  SAlneia  (I) 


Prefaaaer  lym  S.  Baadl* 
Prlt*  Bifinaarlnc  1^. 
Lahlfh  Unlvoralty 
Bathlabaa,  PoanaylTonlJi 


Frofaaaor  B.A,  Belay 
Dapartaant  of  Civil  li«lna*rln< 
Coluabla  Uslvaralty 
Bm  Xork  27,  1,7.  (1) 


Chief,  Duiaau  of  Aaronautlee 
Departaant  of  tha  Navy 
Waahlngttfi  25,  D.C. 

Attni  OB  (1) 

Chlaf,  Bureau  of  Ordnance 
Dapartaant  of  tha  Navy 
Wiahlngton  25,  D.C. 

Attni  Ad3  (1) 


Chlaf,  Bureau  of  Ordnanea 
Dapartaant  of  tha  Navy 
Vaahlngtan  25,  D.C. 


Profaaaor  N.J.  Hoff,  Head 
Division  of  Aeronautical  Eng. 
Stanford  Vnlvaralty 
Stanford,  California 


Preftaaer  V,  Plugga 
Dapartaant  ef  Naehnieal 

luAferd  Onivaralty 
Stanfesd,  CaaifomU  (1) 

Profaitor  X.  Natanyl 
Tha  Ttehnelogleel  Inatlwta 
Northvaitam  Unlvaralty 
Dvanattfi,  Illinois  (1} 


Profaiasr  J.  Xaagner 


Profassor  D.C.  Druelar 
Division  of  KMlaaarlng 
Brawn  UniTaralty  ... 

Provldaaica  IB,  nada  lalaad{Il 

Prefaaaer  J.N,  Ooedlar 
Dapartaant  of  Naohanioal 
XhClnaarthg 
Stanford  Unlvoralty 
StWird,  Callfemln  (1) 

Profaaaor  P.O.  Kadga 
minols  InatlWita  ef  Tachnelegy 
Taehnele<y  CaaiUr  ,  , 

Chicago  16,  nilMls  (1) 


Profaaaor  H.L.  Lmghaar 


Dipt,  of  Aeronaut,  big.,  Apj^.NadilaparWvint  of  ThaoratkaX  aatd 
relyuehnle  Inatlwta  of  Brooklyn  *rpll*l  Maehanlea 
(1)  ?5  LlTlngalen  Svraat  _  UnlvarsUy  of  minols 

L^eoklyn  2,  Nav  Xerk  (1)  Urbana,  minols  (1} 


Profaaser  E.U.  lea 
Division  of  Applied  Mathanatlcs 
Stanford  VnivarelLy 
Stanford,  California  (1} 


Profaaaor  R.D.  Klndlln 

Dapartnaot  of  Civil  &iginaarln|: 

Couunbla  Vnivarelty 

632  W,  125U»  Street 

New  Tork  27,  N.X.  (1) 


Profaaaor  Oaorga  N.  la*  Mr.  N.N.  lancoa 

Dlraotor  of  P^iaarch  Southwaat  Raaaaroh  JbiatituU 

Ranaaalaar  Polytaohnlc  Inallluta  8500  Culabrt  Road 


Profaaaor  WUllan  A*  Naah 

Chlaf,  Buraau  of  Xarda  and  Dooka  Chlaf,  Bura«i  of  Xardt  and  Oaoka  piparUwnt  of  biglnaarlng  Mach, 

(« 

Attni  Coda  D,210  (1) 


Vnlvaralty  of  Plerlda 
OalnaavUIa,  florid* 


Son  Antonio  6,  TMoa 


Dr,  A.  Nadal  Profaaaor  PmI  K,  Na^l 

136  Cherry  Valley  Road  SapartaMt  of  bglnaarlag 

Plttabunh  21,  pMMylvani*  Maehanlea 

(1)  Unlvaralty  ef  CaUfml* 

laibalay,  CallfamU  (1) 

Preftaa*  N.M,  Mawulc  Profaaaor  Aria  Mlli#* 

DtparUaanl  ef  Civil  b«la**rliic  Dapartmnt  of  ClvU  Niglmartag 
Vnlvaralty  of  Ullaals  15  Preapaet  Stmt 

Urbm,  miMls  (1)  Tala  Vnlvanlty 


Chlaf,  Buraau  of  Xarda  and  Oacka  Chlaf,  Buraau  of  Xaida  wd  Oeeka  Chlaf,  Buraau  of  Tarda  and  Dacha 


Profaaaer  I.  RalaaMr 
Dapartaant  of  NMkamil 


Dapartaant  of  tha  Navy  Dapartmnt  of  tha  Navy  Dapartmnt  af  tha  Navy  Pravldanaa  12,  Rhoda  isitns 

Haahington  25,  D.C.  Waahlngton  25,  D.C.  WathlngtoB  25,  0.0. 

Attni  Coda  p,22«  (1)  Attni  Ceda  0-222  (1}  Coda  O-UIO  (1) 

Prof.  Barnard  U,  Shaffar 

Chlaf,  Buraau  of  Xarda  and  ScoIcb  Chlaf,  Buraau  of  Xarda  aid  Dooka  Chlaf,  Buraau  of  Xarda  md  Daaha  2***'’y  ^ 


Provldtnaa  12,  Rhoda  IsltrM  (l)  lhaaaehuaatto  ImiltuU  ef 
Taehnalegy 


PrefMaer  M.  ladsaslv 
Dapartmnt  tf  ItMtwilaa 
RaMoalaar  Palytaatei*  teatitat* 


Attni  Coda  D-U.0  C  (1) 

Comndlng  Offlear  and  Dlraotor 
David  Taylor  Nodal  Baaln 
Washington  7,  D.C. 

Attni  Cod*  61,7  (1) 


(1)  Attni  Coda  D-500 


aiding  Offlear  and  Dlraotor  Ca 


idlng  Offlear  and  Dlraetcr 


New  Xork  Vnlvaralty 
golvarolty  NalghU 
New  Xork  53,  *•*  Xork  (1) 

Profaaaor  A.S.  Valaate* 
Dapartmnt  of  Civil  fitglnaarlng 
Unlvoralty  ef  minols 
Urbans,  minols  (1) 


Cmbrldga  3t,  NMaoehusatta  (1) 

Profoiier  KIL  Stambari  Profaaaer  S.P.  Ilmahsn 

Oivlslcn  of  Applied  IhthaMtks  School  of  InglAaarlng 
Brown  Vnlvaralty  Stanford  Vnivaralty 

Preridanea  12,  Rhoda  Island  (1}  Stanford,  CallfomU 


Prefaaaer  Dam  Taung 
Tala  Uolvaraity 
New  Kavan,  Cennaetlcut 


Or.  John  P.  Brahta 
Dapartmnt  of  ftiglmarlM 
(1)  Vnlvaralty  ef  Califemla 

loa  Angalaa,  California  (1) 


U.S.  Nani  OHnMwa  Laboratory 

Miita  OMi,  Muylaad 

Attni  Taehnieal  Library  (|) 


Attni  Coda  700  (l)  Attni  Coda  7b0 

Oaaaandar  Dlraotor 

U.S.  Naval  Mnama  Uberatory  Katarina  Laboratoiy 
WhiU  Oik,  Maryland  Nav  Xerk  Naval  Shi^ 

Attni  Tach.  Ivaluatlen  Dapt.  (1)  7,  Nav  Tcrt 


)fr.  Martin  Ooland,  Visa  Praaldant  Prof.  J.C.  Camak 


SWthwaat  Raaaareh  Xnatituta 

8500  Cultbra  Read 

San  Antcnlo,  Taxaa  (1) 


Prof.  I.P.  Harrington,  Head 


Dapt*  ef  Civil  foglnaarlng 
Celor^e  State  Vnlvaralty 
Port  CaXIlna,  Celemde  (1) 


Prof.  Bugtna  T.  DiunsUa,  Jr. 


Offlcar-ln^%arga  Dlraotor  Offlear-ln^arga  Dipt,  of  Aarenautleal  Er«lnaarlng  Dapt.  of  Aaronautloal  Eng, 

Naval  Civil  Bnglnaaring  Raaaareh  Naval  Air  bparlmnta  Statlm  yndarvatar  bploakn  r.aaaareb  Dl«  Unlvaraltr  of  Clnolnnatl  frWaim  Un^ralty 

and  Ival.  Ub,, U.S.  Naval  C«atr.  Nav.  Air  Hat.  Cairtar,  NavalXaa  N“f5l5  HavS^Shlgyarxl  Cinoinnati,  Ohio  (1)  Prlneateai,  New  Jaraay  (3 


ItilladalGhla  12,  taiwylvanla 
Attni  MAarUla  Lab.  (1) 


Comandar 

WADD 

VrighV  Pattaraon  APO,  Chlo 
Attni  WVIU03  (1) 

C.O. 

USHMCeU 

KMland  Air  lorea  Baaa 
Alhu(|ttar<]ua,  Nav  Maxloo 
Attni  Coda  20  (1) 


Comandar 

WAUD 

Vrlght  Pattaraon  APB,  Ohio 
Attni  WJRKlO  (1} 


Profaaaor  V.  J.  XaU. 

Lapt.  of  Civil  Wdnawrlng 
Vnlrarsliy  of  miMla 
Urbana,  Illinois  (1) 


CoivMnktr 

WAIO 

Vrlght  PattaramAPB,  Oii*,, 
Attni  WWW  ni) 

Laglslatlr*  Rafaranea  Sarvka 
Library  of  Congrcai 
Waohlngtan  25,  D.C, 

Attni  Dr.  K.  Wink  (1) 


